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Abstract 

In this paper we explore the decay r/^ — )■ r+r^ as a probe for a hght pseudoscalar or 
a light axial vector state. We estimate the standard model branching ratio for this 
decay to be ~ 4 x 10~^. We show that considerably larger branching ratios, up to the 
present experimental limit of ~ 8%, is possible in models with a light pseudoscalar 
or a light axial vector state. As we do not include possible mixing effects between the 
light pseudoscalar and the rj^,, our results should be reliable when the pseudoscalar 
mass is away from the ?7b mass. 
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1 Introduction 



It is widely anticipated that physics beyond the standard model (SM) or new physics 
(NP) will be discovered soon at experiments such as the LHC This NP might contain 
new gauge bosons, additional Higgs bosons beyond the SM Higgs, or new quarks and 
leptons. It is generally believed that these new particles will be heavy with masses 
from the weak scale ~ 100 GeV to a TeV. However, light scalars and vector bosons 
with masses in the GeV range or even lower are not ruled out. For instance, light 
scalar states coming from a primary higgs with non SM decays can be consistent 
with existing experimental constraints [1]. One of the ways to probe these light 
states is to look at decays of particles with masses in the 10 GeV range such as the 
T. Data from the present and future B factories can be used to search for these 
states and/or to put constraints on models that predict such states. 

The pseudoscalar bb bound state in the IS configuration, the rjh, was recently 
observed. Two research groups in BaBar observed it in two different experiments. 
First, it was seen in the decay of T(35') — > 7775 [2] with a signal significance greater 
than 10 standard deviations (a). The rji, was observed in the photon energy spec- 
trum using (109 ± 1) million T{3S) events and the hyperfine T(IS') — 775 mass 
splitting was measured to be 71.4;!;3;f (stat) ± 2.7(syst) MeV from the mass m(r]h) = 
9388.9^2:3 -(stat) ± 2.7(syst) MeV.' Soon after, it was also seen in T{2S) ^rji, 
[3] by another group in BaBar, and the hyperfine mass splitting was determined 
to be 67.4;t6(stat) ± 2.0(syst) MeV from the mass m(r/b) = 9392.9^^8 (stat) ± 
1.9 (syst) MeV . In the past, since the discovery of the T{nS) resonances [4] in 
1977, various experimental environments [5, 6, 7] have been used to seek the ground 
state Tji, but without success. Many theoretical models have attempted to predict the 
mass of r]b. Lattice NRQCD [8, 9] predicts the hyperfine splitting to be Ej^j^ = 61±14 
MeV and correspondingly the mass to be m^^ = 9383(4) (2) MeV which is in agree- 
ment with the experimental results. The calculations of perturbative QCD based 
models [9, 10] predict the hyperfine splitting to be E^^^ = 39 ± ll{th)tl{Sas) 
MeV which is smaller than the measured values. Experiments at BaBar have also 
searched for a low-mass Higgs boson in T(3S') — )■ 7/1", — )■ t~^t~ [11] with data 
sample containing 122 million T(35') events. In the same analysis, constraint on 
the branching ratio for rjb — )■ r+r^ was reported as B7l{rih — )■ r+r^) < 8% at 90% 
confidence level (C.L.). 

In this paper we will be interested in probing light scalar and spin 1 states via rjb 
decays. As the rji, is a pseudoscalar a light pseudoscalar and a spin 1 state with axial 
vector coupling can directly couple to rjb. We will assume the pseudoscalar to couple 
to the mass of the fermion as is usually the case for Higgs coupling to fermions. 
Hence, the rjb which is a bb bound state has advantages over the rjc and rj/r)' mesons 
which are cc and qq{q = u, d, s) bound states, respectively. The 775 is expected to 
be a sensitive probe of a light axial vector state. This follows from the fact that 
the longitudinal polarization of the axial vector, ~ k'^, when k'^ the momentum 
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of the vector boson is much larger than its mass. Consequently, the effective axial 
vector-fermion pair coupling is proportional to the fermion mass for the longitudinal 
polarization. 

In this work we will study the process rjb — > t+t~ mediated by a pseudoscalar 
(A^) or an axial vector (U). In the SM this process can only go through a Z exchange 
at tree level and is highly suppressed with a branching ratio ~ 4x I0~^. There is also 
a higher order contribution to rjb — )■ r^r~ in the SM, via two intermediate photons. 
The branching ratio for this process is also tiny ~ 10~^°. Hence, a measurement of 
BR[rii, — > T~^T~] larger than the SM rate will be a signal of new states. One can also 
probe the states A^{U) in T decays. To search for light A^{U) states in T decays one 
generally considers the decay chains, T A^{U)'y {A^{U) t~^t~) [11]. In other 
words, the A'^{U) is assumed to be produced on-shell. One then looks for a peak 
in the invariant mass of the r pairs. The experimental measurement/constraint 
of BR[T — )■ A'^{U)'~f] X BR[A^{U) — ?■ r+r~] can be converted into a measure- 
ment/constraint on the coupling of the A^{U) to bb, and hence on model parameters, 
if the BR[A^{U) — )■ r+r^] is used as an input [12]. Clearly as niAo^u) > '^r, the 
A^{U) can no longer be produced on-shell and the rate for T — > t'*"t~7 will fall and 
consequently the constraints on the model parameters will be weaker. Note that 
the constraint m^o < 2mB needs to be assumed in the very particular case where 
the CP-even Higgs mass nih < 114 GeV and h 2A^ dominates over h 2mb [1]. 
In general m^o > 2mB is also possible. We will just assume the existence of light 
pseudoscalar and axial vector states close to the rji, mass but they can have masses 
that are greater than or less than 2mfc. 

The r}b has only been seen in the radiative decays T jrjb. Hence, the decay 
rjb — > r+r" has only been studied via the decay T — > t'^t~^. However, the decay 
r)b — >■ T'^T~ can be studied independently from the process T — >■ r'^T~7 as the rjb can 
be produced from various other processes such as two-photon collisions, 77 — )■ r)b 
[6], and in two parton collisions [7, 13], in hadron colliders like the Tevatron and the 
LHC. The process rjb t^t^ has several advantages over T decays in probing A^{U) 
states specially when A^{U) is off-shell which is always the case when itiao^u) > mr- 
First, unlike the rjb which can couple directly to A'^{U), the T can only couple to 
A^{U) in conjunction with another state- usually a photon. Hence, the T couplings 
are second order and therefore it can decay only to the r+r~7 state with a rate much 
smaller than the rate for rjb — )■ r+r". However, the T states are narrower than the 
r)b, which may compensate partially the larger rate for rjb — > t'^t~ relative to T — > 
T~^T~'y in the branching ratio measurements. Secondly, an important distinction 
between T — )■ r^r~7 and rjh — )■ t~^t~ is that the former decay can also proceed as 
a radiative decay in the SM while the latter decay is highly suppressed in the SM 
as indicated above. Adapting the expression used to estimate the SM branching 
ratio for J/ip ^ e+e~7 [14], with the 7 emitted from the final state electrons, to the 
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decay T — >■ T~^r 7, the rate for this decay in the SM is, 



1/3 



2a dE'^ s' 1 - cos^e'^^ 
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Here E'^ represents the 7 energy, 9'^ and the 7 angles, and 9'^ and (l)'^{fl'^) 

the r angles, all in the r+r" cm. frame. is the r velocity and ^^^^ is the angle 
between the r and 7 directions, also in the t~^t~ cm. frame while s' is the r+r~ 
invariant mass squared and s is the T invariant mass squared. The parameter 
A is determined from the experimental data to be (0.88 ± 0.19) [14]. Using the 
branching ratio for T -> r+r" = 2.6 x 10^^ [15] we estimate the branching ratio for 
T ^ r+r-7 = 4.4 X 10-3 with E^ > 100 MeV. 

Naively, the rate for T — > r^r~7 through an off-shell A^, from a 2HDM of 
type II, relative to the SM rate for T — >■ t+t~7 is ~ ^ Tel-^M^''"^^ - Therefore, for 
large tan ^5 ~ 28 the SM and the NP rates may be comparable. However given the 
hadronic uncertainties in estimating the SM and the NP rates for T — )> t~^t~^, it 
will be difficult to distinguish between the NP and the SM contributions. Hence, 
searching for A^{U) with rriAo^u) > '"^t in T ^ t^t'^j will be very difficult because 
of the large SM background. Note that even in e+e~ machines like the S-factories 
where the rji, is produced through the decay T j-qb, the product of branching ratios 
BR[T — )■ 7?7fe] X [?7b — )■ T~^T~] is tiny in the SM because of the highly suppressed 
BR[r]h T+T-] ~ 4 X 10-9. Using the measured BR[T -fr]h\ 10-^ [2, 3] one 

obtains BR[T — )■ 777b] x [rjb — >■ r+r"] ~ 2 x lO"-*^^ which is very difficult to measure. 
In the presence of new physics this product of branching ratios is enhanced and 
can reach < 10~^. Hence the observation of T — )■ 'jt'^t' , with the r pairs coming 
from rjb, at branching ratios much larger than the SM expectations will be signal 
for new light states. In summary, the large SM background in T ^ r+r-7 a 
tiny SM contribution to rjb — > t'^t~ makes the later decay potentially a better probe 
for A^{U) than the former if the decays proceed through the off-shell exchange of 



There are good theoretical motivations for the existence of a light CP-odd 
Higgs boson or an axial vector boson U with masses, m^o and mu respectively, 
in the GeV range or below. There has been interest in the m^o < 2mB region, 
for which a light Higgs, h, with SM-like WW, ZZ and fermionic couplings can 
have mass rrih ~ 100 GeV while still being consistent with LEP data by virtue of 
h — )■ A^A^. This scenario could even explain the 2.3 a excess in the e'^e~ Z + 2b 
channel for ~ 100 GeV [16]. Such a light pseudoscalar Higgs can naturally arise 
in extensions of MSSM with additional singlet scalars and fermions (gauge-singlet 
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supermultiplets) known as Next-to-Minimal Supersymmctric Model (NMSSM) [17]. 
Constraints on models with a light state have been studied recently within a 
2HDM framework with certain assumptions about the couphng and in NMSSM 
[12, 18, 19]. 

Our goal will not be to work in a specific model but we will assume the couplings 
of the to the b quark and the r lepton to be the same as in the 2HDM. We 
will assume this 211DM is part of some extension of the SM. Hence, we will not 
strictly follow the bounds and constraints obtained in some specific extension of 
the SM which includes the 2HDM, but will choose values for the parameters in 
our calculation which are similar to constraints on these parameters in specific NP 
models. The process rjb — )■ t~^t~ will proceed through an off-shell A^ and we will 
consider both m^o < and m^o > m^^. In general, there will be mixing between 
A^ and the rjb and as the pseudoscalar state gets close to the rjb mass the mixing 
between the states will become important [20]. The calculation of this mixing is 
model dependent and while there are estimates of this mixing in simple quark models 
the mixing may be very different in other approaches to the bound state problem 
in QCD. Hence, we will not take into account mixing in our analysis. Therefore, 
our results will be reliable when the mass is away from the rjb mass. We will 
further assume that the A^ is narrow and neglect its width in our calculations. This 
approximation will be good as long as m^o is sufficiently away from the rjb mass. 
When A^ is produced on-shell both mixing and width effects will become important 
and our results will not be reliable. 

There are also models, for example within SUSY with extra gauged U{1), which 
have a light axial vector state [21]. These light states can also mediate the process 
rjb — >■ T+T~. Constraints on these models have been studied [22, 23, 24, 25, 26]. We 
will consider rjb — > t+t~ through the exchange of the axial vector U. To perform 
our calculations we will choose the model discussed in [23, 26] and neglect the width 
of the f/-boson. 

Finally, we note that there are recent dark matter models [27] that also contain 
light scalar (pseudoscalar) and vector (axial vector) states which may be probed via 
rjb — )■ r+r". The HyperCP collaboration has some events for the decay S"*" — >■ pjx'^ji" 
which may be interpreted as evidence for a light pseudoscalar state [28]. 

This paper is organized in the following manner. In section 2 we perform the 
calculations of the decay rjb — )■ t^t~ in the SM and in models with a light pseu- 
doscalar A^ and a fight axial vector U state. In section 3 we present the numerical 
results of the branching ratios for rjb — >■ t'^t" . Finally, in section 4 we present our 
conclusion. 

2 r]b^ T+T- in the SM and NP 

In this section we will study r]b — >■ t+t~ in the SM and in models of NP. The rjb is 
a pseudoscalar and cannot couple to 7 directly. Hence, in the SM, rjb — >■ t'^t~ can 
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only proceed through the exchange of a Z at tree level and we will calculate the 
branching ratio for this process in the SM . This decay can also proceed at higher 
order in the SM through intermediate two photon states. 

In the presence of NP rjb — > r+r" can proceed through the exchange of a hght 
pseudoscalar or a light spin 1 boson with axial vector coupling. We will consider 
these two NP scenarios in this section. The various tree level contribution to the 
rjf, — )■ r"^r~ in the SM and NP are shown in Fig. 1 and Fig. 2, respectively. 




Figure 1: Various processes contributing to rjb t^t in the SM. 




A" 





Figure 2: Various processes contributing to rji, — > r+r in NP. 

We begin with rib — )■ t'^t~ in the SM. We show, in Fig. 1, the decay process 
77b — > T^T~ via the Z-boson exchange and through the two photon intermediate 
states. The decay rate for the tree level Z exchange process can be obtained as, 

,2 



r^(r7b ^ r+r-) 



167r cos^ 6w 

where 9w denotes the Weinberg angle, = ^1 — ^^^^ is the velocity of the 
lepton in the r]b rest frame and 

1 



(3) 



K -M|)2 + M|r|- 



The decay constant /^^ in Eq. 3 is defined as [29] , 



(4) 



(5) 



The process r]^ t~^t~ can also go via two photon intermediate states as shown 
in diagram Fig. 1. This diagram is dominated by the imaginary part [30] which we 
can estimate using unitarity [31] to obtain. 



, (1 + /3. 
In — 



m 



Vb 



T[rjb 77], 



(6) 



where a is the electromagnetic fine structure constant. One can calculate V\r]}, 77] 
as, 



81m? 



r[r/.^77] = o.:;r" . (7) 



where we have used the heavy quark limit for the h quark. Since the 27 exchange 
contribution is mostly imaginary relative to the Z exchange contribution therefore 
to a good approximation the total width Tt\j]h — >■ r+r"] is , 

Vt\r]^^T^T-\ « r^[r76^T+T-]+r2^[r76^T+T-]. (8) 

We now turn to NP models and begin with the 2HDM. The couplings of the down- 
type quarks D and charged leptons I with in the generic 2HDM model are given 
by [32] 

CPJ = '-^{D]VC^^,D + lMf^'^d)A\ (9) 

where is a model-dependent parameter, M^*"^ = {md,mc,mi,) and Ivlf""^ = 
(me,mfj_,mr) are the diagonal mass matrices of D and £, respectively. We will 
consider F^o > 1 in our analysis. In the case of 2HDM type (II) F40 = tan(3 while 
in 2HDM type (I) F40 = - cot/3. 

Ill Fig. 2(a) we show the decay process rjb — )■ r'^r" via the exchange of the 
CP-odd Higgs scalar A^. The decay rate for this process can be obtained as, 

r^\m^r-r-) = ^^f^^M', (10) 

where the coefficient a^o depends on the mass m^o as, 

la^oP = " (11) 

We have assumed that the decay width F^o for the is negligible. In Eq. 10, we 
have used, 

(0| 6(0)75^0) |r?.(g)) = (12) 

where /^^^ has been defined in Eq. 5. 

Finally, we move to NP models that contain a light spin 1 boson with axial vector 
couplings. In Fig. 2(b) we show the decay process rjb — > t'^t" via the exchange of 
the light neutral gauge boson U. We write down a model independent Lagrangian 
for the U -boson but we assume the structure of the Lagrangian to be similar to the 
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one discussed in Ref. [23, 24, 25]. We take the U couplings to the down-type quarks 
and charged leptons to be given by 

j^D,i ^ f^'\DYl,D + lr^d)U^, (13) 

with the axial coupling 

f^^' = 2-lGlmuFu, (14) 

where mu denotes the mass of ^7-boson and Fu denotes a model-dependent param- 
eter. In the specific model [23, 24, 25], Fu = cos C tan /3. 

Again, we will be interested m. Fu > 1. The decay rate for rfb — > t"^t~ can be 
obtained as 



r^(r7.^r+r-) = ^^^^^-^^ ^^(^2^ _ ^j2^4|^^|2^ ^^^^ 



where 



Eq. 16 can be expanded as. 



«t^r = 7-^— i2T¥(i-^' +•••), (17) 



{ml - ml) 



if 2^ — 7i — ^2 y '1 \ ^ 1 



Neglecting x, Eq. 15 reduces to 



V\n,^T^T-) = (18) 

iDTT 

Thus, Eq. 18 shows that the decay width for r\\, t^t~ does not depend on mu 
in the approximation of neglecting the width of the [/-boson. This result is easy 
to understand. If one increases the mass of the U then the matrix element for 
T^T^ is suppressed due to propagator effects. However, the coupling, which 
is proportional to mu, increases to compensate for this suppression. The fact that 
the width for 775 — > t^t~ is independent of mu only holds because the r]b is a 
pseudoscalar. 

The result of Eq. 18 does not make sense as mu gets sufficiently large as the 
couplings in Eq. 14 becomes non-perturbative. Requiring the couplings to be < 1 
one gets the constraints mu < Hence for Fu ~ 50 one can get mu to be in 

the GeV range. 
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It is interesting to note that in the up sector the behavior for the decay width is 
different. The coupling of the vector boson to the up type quark, U, is given by 



(19) 



with the axial coupling of the up-type quarks 



(20) 



In the model of Ref. [23, 24, 25], = cos (cot/?. 

For instance, the branching ratio BTZ{r)c — > fJ''^IJ'~) does not depend on mu or on 
tanf3 and is given as. 



where = y 1 — y^^j and /^^ is the ijc decay constant. We can see from Eq. 21 

that the branching ratio BTZ{ric — )■ /i'''/i~) is much smaller than BTZ{r]i, — >■ r'^r~) if 
tanP > 1 because of the absence of the factor tan^ (3 in the rate for rjc — >■ A*"*"/^" ■ 

3 Numerical Analysis 

In this section we present our numerical results. We take the average rib{lS) mass 
to be rrirj, = 9390.8 ± 3.2 MeV [3], the decay constant = (705 ± 27) MeV [33] 
and the width to be T^^ 10 MeV [34]. 

In the SM, at tree level, rjh — )■ r+r" goes through the exchange of a Z-boson and 
we obtain a tiny branching ratio BTZ^{7]b — )■ r+r~) = 3.8 x 10"^. In our calculation 
we have used Tz = 2.4952 ± 0.0023 GeV [15]. For the two photon contribution to 
rjb — >■ T~^T~, we obtain, using Eq. 6 and Eq. 7, BTZ^^lrib — >■ r+r^] > 4.6 x 10^^° for 
ruf, = 4.8 GeV. Using Eq. 8 the total branching ratio for rjb — )■ r+r^ is ~ 4.3 x 10"^. 

In Fig. 3, we plot the logarithm of the branching ratio for 1]^ — )> r"^r" mediated 
by the pseudoscalar in a generic 2HDM model. The branching ratio, BTZ^° , is 
plotted for various values of the ^4° mass, which we take from 0.1 to 20 GeV, and 
for various values of F^o. As the mass of the ^4° approaches the mass of the rjb 
the branching ratio increases and blows up at rrij^o = m^^^. This behavior clearly 
does not represent the physical situation because in this region the width of the 

and mixing effects of the A^ with rjb become important and regularize the A^ 
contribution. We observe in Fig. 3 that the branching ratio ~ F^o is very sensitive to 
Fj^o. The branching ratio is relatively less sensitive to the mass m^. We see from the 
plots in Fig. 3 that the branching ratio for rjb — ?■ r+r^, through the exchange, can 
be considerably larger than the SM branching ratios and can vary from ~ 10~^ to 
the experimental bound of 8 % for F^o — 40. Since we have neglected the width and 



r^(7;e^/xV) 



167r 



■j3r COS^ 



(21) 
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U , , , , , , , , , , , , , , , , , , _|qU , , , , , , , , , , , , , , , , , , , ^ 

2 4 6 8 10 12 14 16 18 20 

mj(i[GcV] m^(i[GcV] 

Figure 3: The logarithm of BlZ^^\'qij t^t^) as a function of rriy^o for different 
values of F40 and m^o G [0.1,20] GeV. 

mixing effects our predictions are no longer reliable as the mass of the approaches 
the mass of the r]b. The mixing effects are model dependent and as an example, for 
the model for mixing employed in Ref. [20], the effects of mixing are important in 
the m^o mass range of 9.4 — 10.5 GeV. We see from Fig. 3 that even outside this 
range the branching ratio for rjb — ?■ r+r~ can be significant and we expect the same 
to be true also in the mass range where mixing effects are important. 



-2 - 




_1qU|^ : : : , : : : : , : : : : , : : : : U 

10 20 30 40 

Fu 

Figure 4: The logarithm of BTZ^ {7]^ — > r+r^) as a function of Fjj. 

As discussed in the previous section, the branching ratio for the decay BTZ^ {r]h — > 
T^T^) is independent of the mass of the gauge boson U in the approximation of 
neglecting the width of the f/-boson. We next plot in Fig. 4 the logarithm of the 
branching ratio for r/;, — )■ r+r" versus Fu. Working in a specific model [23, 24, 25] 
Fjj = cos C tan /3. We plot the branching ratio versus the invisibility factor cos^ for 
different values of tan /3 in Fig. 5. Again we observe that the branching ratio can 
vary over a wide range and can be much larger than the SM prediction. 
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Figure 5: The logarithm of BTZ {rjh — t- r+r ) as a function of cos^ for different 
values of tan (3 and cos ( G [0, 1]. 

4 Conclusion 

In this paper we explored the decay rjb — )■ r+r^ as a probe for a light pseudoscalar 
or a light axial vector state. We estimated the SM branching ratios for rib — )■ t~^t~ 
via the Z exchange and the two photon intermediate state and found it to be very 
small ~ 4 X 10~^. We then considered the decay process i]b — )■ r+r^ mediated via 
the pseudoscalar Higgs boson in a 2HDM type NP model. We found that the 
branching ratio for rjiy — )■ t~^t~ can be substantially larger than the SM prediction 
and can reach the experimental bound of 8 %. Working in a specific model containing 
a light axial vector state, U, a similar result was obtained for the branching ratio 
of rjb —7- r+r~. We also obtained an interesting result that the BTZ^ {rjb — t- r+r^) 
is independent of the mass of ?7-boson if the width of the U is neglected. This 
result followed from the fact that the axial ?7-boson couplings to fermions were 
proportional to the mass and the fact that rjb is a pseudoscalar. A constraint on 
the [/-boson mass could be obtained by requiring its coupling to fermions to be < 1. 
In light of the results obtained in the paper an experimental measurement of the 
branching ratio for rjb — )■ t~^t~ is strongly desirable as this measurement might reveal 
the presence of light, ~ GeV, pseudoscalar or axial vector states. The experimental 
measurements of rjb — t- r+r^ may be feasible at planned high luminosity B factories 
and at hadron colliders such as the Tevatron and the LHC, specially if the branching 
ratios are much larger than the SM rate. 
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